Interacting Hofstadter spectrum of atoms in an artificial gauge field 
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Motivated by experimental advances in the synthesis of gauge potentials for ultracold atoms, we 
consider the superfluid phase of interacting bosons on a square lattice in the presence of a magnetic 
field. We show that superfluid order implies spatial symmetry breaking, and predict clear signatures 
of many-body effects in time-of-flight measurements. By developing a Bogoliubov expansion based 
on the exact Hofstadter spectrum, we find the dispersion of the quasiparticle modes within the 
superfluid phase, and describe the consequences for Bragg spectroscopy measurements. The theory 
also provides an estimate of the critical interaction strength at the transition to the Mott insulator 
phase. 



The spectrum of a particle on a tight-binding lattice in 
the presence of a magnetic field [l|-|3| is a problem that 
is simple to state but has surprisingly rich phenomena. 
In the infinite system, it is sensitively dependent on the 
precise value of a, the magnetic flux per plaquette of 
the lattice (measured in units of the flux quantum). For 
rational a = p/q (p and q coprime), the spectrum splits 
into q bands, and each state is g-fold degenerate. When 
the density of states as a function of energy is plotted 
against a, the resulting 'Hofstadter butterfly' [3[ has a 
fractal structure. 

Recent work using cold atomic gases has raised the 
possibility of observing the Hofstadter spectrum directly 
in experiment 043 ■ Two distinct approaches to produc- 
ing effective magnetic fields for neutral atoms have been 
proposed and implemented. In a rotating optical lattice 
[8Hl0l|. the Coriolis force simulates the magnetic Lorentz 
force, while external lasers applied to a stationary lat- 
tice can imprint motion-dependent phases to synthesize 
a gauge potential [III. \12\. 

Previous theoretical work on this system [lj| has con- 
sidered vortex pinning for both shallow [l]| and deep Q 
lattice potentials, and the transition from the superfluid 
to a Mott insulator 14 - 1|| . In this work, we introduce an 
alternative theoretical approach to the system, using Bo- 
goliubov theory [l!| [2(| to provide a controlled expansion 
for the superfluid phase, in terms of both thermal and 
quantum fluctuations. We consider the case with both 
a and mean density p (particles per site) of order unity, 
which is of direct relevance to experiment. To the best 
of our knowledge, this is the first theoretical study of the 
superfluid phase in this parameter regime using the Bosc- 
Hubbard model. Starting from this microscopic descrip- 
tion, and including the exact Hofstadter spectrum, we 
calculate the condensate configuration based on a mini- 
mization of the on-site interactions, without resorting to 
a phenomcnological Ginzburg-Landau functional. This 
aspect is similar to earlier work on frustrated Josephson 
junction arrays 21-22 



| , and related models of bosons and 
fermions in optical lattices 24 , 2f| . 

Besides describing the real-space condensate configu- 
ration, which demonstrates the intricate interplay of the 



magnetic vortices and the external lattice potential, our 
approach predicts features that should be directly observ- 
able in time-of-flight images. In particular, we show that 
the superfluid necessarily breaks translation symmetry, 
resulting in extra peaks in time-of-flight images, clearly 
distinguishing this case from that with zero flux (2(| and 
also with noninteracting bosons . 

We then go beyond mean-field theory to calculate the 
spectrum of Bogoliubov quasiparticlcs, which inherits the 
complex structure of the Hofstadter spectrum, while also 
exhibiting the characteristic features of the Bose-Einstein 
condensate, including a linearly dispersing phonon mode. 
This spectrum is accessible in experiments using Bragg 
or lattice-modulation spectroscopy, and we make predic- 
tions for the appropriate response functions. 

The approach that we describe is valid deep within the 
superfluid phase, in contrast to previous theory 
and is directly relevant to experiments. Within our the- 
ory, we calculate the depletion of the condensate, due to 
both thermal and quantum fluctuations, providing lim- 
its on the validity of the approach and also approximate 
boundaries for the superfluid phase. The range of appli- 
cability can be extended by taking into account higher- 
order terms in the Bogoliubov expansion, which give in- 
teractions between the quasiparticles. 

We begin with the Hamiltonian for a single-band Bose- 
Hubbard model on a two-dimensional square lattice, in 
the presence of a uniform magnetic flux of a per plaquette 
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where bj and rij are the annihilation and number oper- 
ators on site j with coordinates (xj,yj); t, U, and fi are 
the hopping strength, Hubbard interaction, and chemical 
potential; and u = e 2ma . We use the Landau gauge, with 
the vector potential vanishing on links in the x direction, 
which is both theoretically convenient and best suited to 
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experiments using an optically-induced gauge potential 
fill Il2j . For rotating lattices the symmetric gauge 
is more natural, requiring a straightforward gauge trans- 
formation to be applied. 

For simplicity, we will treat an infinite spatially- 
uniform system and assume rational a throughout. The 
main effect of the finite trap size is to wash out the small- 
scale fractal structure of the Hofstadter butterfly 0,0): so 
our results should be valid for either rational or irrational 
a. We assume that other effects of the trapping potential 
can be incorporated using a local-density approximation. 

Assuming rational a = p/q, one has uj q = 1, so the 
magnetic unit cell spans q sites in the x direction, and 
the Brillouin zone (BZ) is correspondingly reduced. It is 
helpful to introduce the magnetic translation group 
with primitive translation operators T x , v obeying T x bj - 
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and Tybj 



former expression, resulting from using the dynamical 
rather than canonical momentum, gives [7 x ^ y ,'H] = 0. 
In momentum space, T x shifts by Y = 2nay, while the 
hopping term in % mixes states separated by X = 2nax. 

This leads to a noninteracting spectrum Q consisting 
of q bands (labeled by 7), within which each state is 
q-fold degenerate. The single-particle states ip in {k) are 
given by linear combinations of momenta fc + nX for 
n = 0, . . . , q — 1. We define annihilation operators 
for crystal momentum k -\-lY , where I = 0, . . . , q— 1 and 
fe is the momentum referred to the doubly-reduced BZ 
23 , — ^ < k x ,k y < ^ . The energy associated with these 
single particle states is e 7 (fc), independent off. (The case 
a = 5 is shown by the dashed line in Figure [1] below.) 

Written using the operators aki-y, the kinetic term in % 
is diagonal, while the interaction term can be expressed 
using coefficients u depending on four sets of fc, £, and 



7. These coefficients are constrained by symmetries [28 1 
such as translation, which enforces conservation of crystal 
momentum. 

Within the superfluid phase, the U(l) phase symme- 
try is spontaneously broken and the boson annihilation 
operator has a nonzero expectation value Following 
Bogoliubov [2^ , we use the ansatz 



a-ki-t = A tl (2tt) S (fc) + a k e~/ 



(2) 



where the c-numbers give the condensate order pa- 
rameter and the operators dke-y describe fluctuations out 
of the condensate. The condensate can contain contri- 
butions from all bands 7 and from crystal momenta £Y 
with nonzero £, and so also breaks symmetry under T v 
(but not T3)- Note that, at least within our mean-field 
theory, translational symmetry is simultaneously broken 
at the condensation transition Hi} . 

Our approach is to expand the Hamiltonian in pow- 
ers of the operators aki-y, which in physical terms is an 
expansion in fluctuations, both thermal and quantum. 
The lowest order in this expansion contains no opera- 
tors and leads to mean-field theory, while the terms with 



a single operator cancel. The quadratic terms describe 
Bogoliubov excitations above the condensate, and higher 
orders give interactions between these modes. 
The mean-field energy density is 



bj+yTy The phase factor in the (bj) — Y^in 
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where u is the interaction coefficient u with all mo- 
menta fci...4 = (and a factor of volume removed). The 
quantity ho should be minimized with respect to Ag~ t 
to give the mean-field condensate wavefunction in mo- 
mentum space. The corresponding real-space quantity, 
XL ip-y n (0)Ae-y, is a function of 



n Xj +£ yj 



and Dj modulo q; the ansatz of Eq. ^ implies a real- 
space unit cell of q x q sites. (The density pattern min- 
imizing the mean-field energy for a = 5, which consists 
of diagonal stripes, is shown in the right inset of Figure Q] 
below.) The minimization of h$ is therefore equivalent to 
minimization with respect to the real-space condensate 

14[ is more 



configuration. The latter perspective [1 
appropriate in the continuum, whereas here the lattice 
provides a strong pinning potential that simplifies the 
momentum-space approach. 

The q x q real-space unit cell implies a reduced BZ, with 
observable consequences in time-of-flight images, where 
Bragg peaks are expected at points corresponding to mo- 
menta fc equal to reciprocal lattice vectors [2(|. (Note 
that the natural gauge for timc-of-flight measurements 
depends on the experimental procedure. We use the Lan- 
dau gauge appropriate to current experiments with artifi- 
cial gauge potentials ll|,|l2j.) The peak for fc = nX+lY 

2 , multiplied by an envelope 



has intensity | A ej tp yn 
due to the momentum-space Wannier wavefunction. We 
find that on-site interactions generally favor condensate 
configurations with intensity (but not phase) indepen- 
dent of £, implying that the only I dependence of the 
peak intensity results from the Wannier envelope. Sig- 
nificantly, the peaks with £ ^ result from spontaneous 
translational symmetry breaking and are not expected 
in the noninteracting case 0, so time-of-flight measure- 
ments give a clear signature of many-body effects. 

Like the single-particle states, the condensate config- 
urations are at least g-fold degenerate, corresponding to 
different symmetry-related spatial orderings. This allows 
for the possibility of real-space domain formation, which 
would not affect time-of-flight images and would likely 
require more sophisticated in situ probes to confirm. 

Minimization of the mean-field energy ho causes the 
terms in % linear in the operators aki-y to vanish. 
Combining the creation and annihilation operators into 



and 'flattening' this into a 2q 



component vector for each fc, the quadratic terms can 
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be written as 



H V) = Z 



fees 



d 2 fc 



alM(k)a k +n c 
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where H c is a constant. Similarly to /io, the matrix M(fc) 
has contributions from both the kinetic and interaction 
energy; the latter are self-energy terms for the quasipar- 
ticles due to scattering with bosons in the condensate. 
They include 'anomalous' processes in which a pair of 
condensed particles scatter from each other into an ex- 
cited state and the reverse process where they return to 
the condensate, resulting in quadratic terms that do not 
conserve the number of a*^ 7 quanta 1|J 20| - 

In the limit k — > 0, M(0) is the Hessian matrix for the 
mean-field energy ho (with respect to variations of 
and its conjugate) and so is a nonncgativc-definitc ma- 
trix, with a single vanishing eigenvalue corresponding to 
the U(l) symmetry of ho- For nonzero fc, all eigenvalues 
of M(fc) are strictly positive. 

To find the quasiparticle modes, one must diagonalizc 
the quadratic form in Eq. (0J using a Bogoliubov trans- 
formation, which for bosons amounts to finding the sym- 
plectic transformation that diagonalizes M(fe) [29|. The 
quadratic part of the Hamiltonian then becomes 



H (2) = 



ke<s 



d 2 k 
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where dk( is the annihilation operator for the Bogoli- 
ubov mode labeled by ( = l,...,g 2 , with energy 
While the Bogoliubov quasiparticles have definite mo- 
mentum referred to 25, they cannot be labeled by £, be- 
cause the condensate configuration breaks the symmetry 
under translation T y . 

The constant H' c in Eq. ((5|) gives the zero-point cor- 
rection to ho which, along with the contribution from 
thermally excited quasiparticles, can in principle change 
the global minimum of the free energy and hence the 
most stable condensate configuration. 

Using the fact that M(fc) is positive definite for k ^ 0, 
one can show that the mode energies are all real and 
positive [29[ . At k = 0, one of the eigenvalues vanishes, 
giving the Goldstone mode resulting from the broken 
phase rotation symmetry. This mode has linear disper- 
sion near fe = 0, with phonon velocity independent of the 
direction of fc, while the remaining q 2 — 1 modes have a 
nonzero gap and quadratic dispersion. Figure [1] shows 
the dispersions £kc for a = 5, found by minimizing ho 
and diagonalizing M(fe) numerically. For all a, there are 
q 2 modes, including one Goldstone mode; here only the 
lowest 6 are shown, for clarity. 

The Bogoliubov ansatz, Eq. ([2]), and expansion in pow- 
ers of operators is in principle exact, with the higher- 
order terms leading to interactions between the quasi- 
particles. The present approximation, truncating the se- 
ries at quadratic order, is valid provided that the system 




FIG. 1: (Color online) Quasiparticle dispersion (solid lines) 
and noninteracting single-particle dispersion (dashed), both 
in units of hopping t, for a = | = |. The dispersions are 
plotted along a path in the reduced Brillouin zone 23 shown 
in the left inset. In the interacting case, U = 2t, the mean 
density is p = 1, and the real-space density pattern is as 
shown in the right inset (one unit cell of q x q sites; black 
squares are sites with higher density than gray squares). In 
both cases, there are q 2 = 9 modes, of which only the lowest 
6 are shown. For U = 0, the modes are q-fold degenerate and 
have been shifted vertically by an arbitrary choice of chemical 
potential. 



can be treated as a low-density gas of quasiparticles, or 
equivalently, that the system is deep within the super- 
fluid, with only weak thermal and quantum fluctuations. 

This criterion can be quantified by calculating the con- 
densate depletion, the quasiparticle contribution to the 
total particle density. The leading contribution is the 
mean- field term due to the condensate, J"^ |A£ 7 | 2 , while 
the first-order correction includes both thermally excited 
quasiparticles and zero-point quantum fluctuations, inte- 
grated over 23. We define the depletion as the ratio of 
this correction to the mean- field result. 

For nonzero temperature T, the momentum inte- 
gral in fact diverges logarithmically, an instance of the 
Mcrmin-Wagncr-Hohcnbcrg theorem |3l|, which for- 
bids breaking of a continuous symmetry in two dimen- 
sions for T > 0. In the infinite system, there is therefore 
no true condensate, and so the 'depletion' is complete. 
In a real system, the small-momentum divergence is cut 
off at a scale ko, given either by the finite size i? e ff of 
the system in the two-dimensional plane (ko ~ R^s) or 
by a small hopping matrix element t± in the transverse 
direction (ko 



\/2m*t_L, where m* is the effective mass 
in the lowest band). 

Figure [2] shows the depletion for a = | , as a function 
of density, interaction strength, and (in the inset) tem- 
perature. It is small deep within the superfluid phase 
and increases to roughly 25% for the largest values of U 
and T shown. Neglecting cubic and quartic terms within 
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tadtcr spectrum of noninteracting particles. Our theory 
describes the spatial symmetry breaking of the conden- 
sate wavefunction and allows the quasiparticle spectrum 
to be found. We predict clear signatures for time-of-flight 
images and Bragg spectroscopy, which we expect to be 
observable in experiments with ultracold atoms in the 
near future. 

We thank Ian Spielman and Trey Porto for helpful dis- 
cussions. This work is supported by JQI-NSF-PFC and 
ARO-DARPA-OLE. 



FIG. 2: (Color online) Condensate depletion for a — | as 
a function of interaction U/t (main figure) and temperature 
T/t (inset), where t is the hopping strength. In the main 
figure, T = and the densities are p = 1 (top curve), 2 
(middle), and 4 (bottom); in the inset, p — 1 and U/t = 2. 
The depletion is smallest, and hence the approximation best, 
deep in the superfluid phase, with weak interactions, high 
density, and low temperature. For T > 0, the logarithmic 
divergence of the depletion integral is removed with small- 
momentum cutoffs of fco = 0.1 (solid line) and 0.02 (dashed 
line) lattice units, corresponding to an effective system radius 
of 10-50 sites. 



the Bogoliubov theory relies on the assumption of small 
depletion, and so the conclusions presented here are only 
qualitatively applicable for larger values of U and T. 

The depletion calculation also provides a rough esti- 
mate for the boundary of the superfluid phase, at the 
point where the depletion reaches 100%, although the 
approximation of independent quasiparticles is probably 
not valid at this point. For a = i, p = 1 and T = 0, 
this gives an estimate of (t/U) c = 0.08, in reasonable 
agreement with the value of (t/U) c = 0.063 (at the tip of 
the p = 1 Mott lobe) found using the Gutzwiller ansatz 
171 ]. It should be noted that the latter approach, which 
neglects fluctuations within the Mott insulator, generally 
underestimates (t/U) c [32j |. 

Techniques to measure the qu asip article spectrum in 
ultracold atomic experiments [33j,l34J include Bragg spec- 
troscopy 



36j , where a two-photon Raman transition 
is used to measure the dynamic structure factor, a four- 
point correlation function of the boson operators bk- In 
the presence of a condensate, the dominant contribution 
factorizes into the product of the condensate density and 
the two-point correlation function of the quasiparticle op- 
erators dk( ■ The structure factor at frequency oj is there- 
fore given by a delta function 5{uj ± £fc^) at the energy 
of each quasiparticle mode £, allowing the quasiparticle 
spectrum to be measured directly. 

We have presented a theory describing the superfluid 
phase of bosons in a magnetic field on a square lattice, 
and showed how many-body effects modify the Hofs- 
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